Fermi-Luttinger liquid: spectral function of interacting one-dimensional fermions 
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We evaluate the spectral function of interacting fermions in one dimension. Contrary to the 
Tomonaga-Luttinger model, our treatment accounts for the nonlinearity of the free fermion spec- 
trum. In a striking departure from the Luttinger liquid theory, the spectrum nonlinearity restores 
the main feature of the Fermi liquid: a Lorentzian peak in the spectral function on the particle 
mass-shell. At the same time, the spectral function displays a power-law singularity on the hole 
mass-shell, similar to that in the Luttinger liquid. 
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I. INTRODUCTION 

One-dimensional problems play a special role in quan- 
tum many-body theory. In many cases, the reduced di- 
mensionality affords one a deeper understanding of the 
role of interactions in a many-body system. Recent 
progress in experimental techniques has also contributed 
to the increased attention paid to a variety of one- 
dimensional (ID) fermionic and bosonic systems. Ex- 
amples include edge modes of the quantum Hall liquid^, 
carbon nanotubes^, cleaved edge semiconductor wires'^, 
antiferromagnetic spin chains'*, and cold atoms in ID op- 
tical traps^. These developments catalyze the interest to 
the fundamental theory of interacting ID quantum liq- 
uids^^i^ 

The traditional framework for discussing ID systems 
is provided by the exactly solvable Tomonaga-Luttinger 
(TL) model^^"*^. The crucial simplification that makes 
the model solvable is the linearization of the fermionic 
dispersion relation. It was understood early on-'^^ that a 
model with a linear spectrum is an idealized one. The 
goal of this paper is to elucidate the influence of the dis- 
persion nonlinearity on the spectral function of ID Fermi 
systems. 

In the absence of interactions, the spectral function 
is given by Ap{e) = S{e — £,p), where is the single- 
particle energy measured relative to the Fermi level (^p = 

/2m — ep for Galilean-invariant systems). According 
to the Fermi liquid thcory^^ , weak repulsive interactions 
merely broaden the peak in Ap{e) to a Lorentzian, 
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Im ■ 



1 



^p + i/2Tp ' 



(1) 



The broadening originates in the finite decay rate l/2rp 
of the Fermi liquid quasiparticles. The rate can be esti- 
mated with the help of the Golden Rule, 



(2) 



where v is D-dimensional density of states and V is 
the characteristic strength of the short-range repulsive 



interaction. Indeed, quasiparticle relaxation occurs via 
real transitions resulting in the excitation of particle-hole 
pairs. In dimensions D > I already a single pair produces 
a finite decay rate. The amplitude of such process is pro- 
portional to V, hence the dimensionless factor (i^V)'^ in 
Eq. (2). The factor accounts for the corresponding 
phase space volume, i.e. the number of possibilities a 
pair can be excited while obeying the energy and mo- 
mentum conservation laws. 

In ID the situation differs dramatically. Indeed, the 
Golden Rule result Eq. (2) is identically zero for = 1. 
Moreover, it can be shown^'' that in the framework of the 
TL model^^"^^ , the self-energy vanishes on the mass-shell 
in all orders of the perturbation theory. In fact, the TL 
Green function^''' ^® assumes manifestly non- Fermi liquid 
form. In the vicinity of the mass-shell, |e — ^fe| <C \£,k\, it 
reads 
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(3) 



where k = p — pp and ^/j = iwfc is the dispersion re- 
lation linearized near the right (left) Fermi point ±pp. 
The corresponding spectral function of, say, the right- 
moving excitations then takes the form (hereinafter we 
concentrate on the limit of zero temperature^^) 
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(4) 



The exponent 70 in Eqs. (3) and (4) characterizes the 
interaction strength and in the lowest order is given by 
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y{Vo-V2p,), 



(5) 



where v = (27rw)~"'^ is the density of states in ID and 
Vfc is the Fourier transform of a short-range interaction 
potential. 

Unlike the Fermi-liquid Lorentzian Eq. (1), the spec- 
tral function Eq. (4) exhibits a characteristic threshold 
behavior. The edges of the spectral support coincide with 
the single-particle energies ^fc, and the spectral function 
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FIG. 1: (a) Support of the spectral function in (p, e)-plane. 
The parabola represents the mass-shell (e = ^p). For e < 0, 
the mass-shell coincides with the edge of the spectral sup- 
port. For e > 0, however, the mass-shell (dashed lines) falls 
into a broad spectral continuum. In this region, the spec- 
tral function has a peak of a finite width, which is indicated 
schematically by dark gray, (b) The low-energy region near 
the right Fermi point p = pp, with k = p — pp- 



displays a power-law edge singularity on the mass-shell 
e — > ^fc. This singularity is a hallmark of the Luttinger 
liquid^*^ behavior. The particle-hole symmetry of the TL 



model implies that Afc(e) = e), hence the behavior 

of Ak{e) in the particle region of the spectrum e > is 
identical to that in the hole region e < 0. 

The exact solvability of the TL model relies crucially 
on the assumption of strictly linear dispersion relation. 
The purpose of this paper is to examine the effects of 
the dispersion nonlinearity on the spectral function of 
ID spinless fermions. Specifically, we consider a non- 
linear dispersion relation with a positive curvature, and 
approximate the single-particle spectrum in the vicinity 
of the right (R) and left (L) Fermi points by 



± vk - 



2m 



PTPF- 



(6) 



The presence of a finite mass m breaks the particle-hole 
symmetry of the TL model and affects the spectral func- 
tion in the particle and hole regions of the spectrum in 
manifestly different ways. 

The effect of the dispersion nonlinearity on the par- 
ticle region e > is the most dramatic. Rather than 
being the edge of the spectral support, the mass-shell 
e = now falls within a broader continuum, see Fig. 1. 
Consequently, the edge singularity [cf. Eq. (4)] on the 
mass-shell disappears and gets replaced by a peak of a 
finite height. Both the shape and the width of the peak 
appear to be rather different from those in _D > 1, see 
Eqs. (1) and (2). For e > and |e - < k'^/2m we 
found 
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The fc-dependent exponent is given by 
where 

f^k^HVo-V,)^ (9) 

is the exponent governing the power-law divergence of the 
dynamic structure factor*. For short-range interactions 
Hk ^ when A; — > 0, while fJ,2pF = 270, hence fc ^ 
limit of Eq. (8) agrees with Eq. (5). 

The main feature of Eq. (7) is the appearance of a 
finite quasiparticle decay rate, l/2Tfc. To the lowest non- 
vanishing order in the interaction strength it is given by 

(10) 



A finite decay rate emerges only in the fourth order in 
the interaction strength^^. This is because the minimal 
relaxation process involves excitation of two particle- hole 
pairs on the opposite branches of the spectrum. The 
factor (^fc^)"* in Eq. (10) accounts for the corresponding 
phase space volume. Note also that I/t^ oc and 
vanishes in the limit m — > 00 taken at a fixed v, which 
corresponds to the TL model. 



Despite its small value compared to that in higher di- 
mensions, the very emergence of a finite quasiparticle 
relaxation rate in a ID system is a matter of funda- 
mental significance. One may even wonder if it puts 
ID Fermi systems back in the realm of the conven- 
tional Fermi liquid theory. Indeed, in a broad region 
k ~ Cfe^l ^ ('''fc7fe)~^ around the mass-shell, the spectral 
function Eq. (7) is essentially a Lorentzian with the width 
1/t/c < (Tfc7fe)"^ < , see Fig. 2. Restoring aU the fac- 
tors, we find for the immediate vicinity of the particle 
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mass-shell, 



< 



A fA-T" fji^V" 



l/2Tfc 



Furthermore, for not too small momenta, such that 



(11) 



ln(fc/pF) 



> 



1 



(12) 



(this conditfon allows for fc pp at ^ 1) the 
Lorentzian peak (11) carries most of the particle's spec- 
tral weight, see Section VII. This is the hallmark of the 
Fermi liquid^^. The Luttinger liquid behavior is found 
only at smaller momenta and sufficiently far away from 
the mass-shell, at \e - > {Tkjl)~^- 

The dispersion nonlinearity has a strong effect on the 
spectral function despite being irrelevant in the renormal- 
ization group (RG) sense: it does not affect, for example, 
the power-law asymptote of the local tunneling density 
of states at low energies. However, momentum-resolved 

measurements'^'^ and numerical simulations^ may reveal nomenon'^'^ . A detailed theory of the threshold behavior 



along the entire hole mass-shell line, see Fig. 1. The 
quadratic dependence of the exponent on the interaction 
strength originates in the orthogonality catastrophe phe- 
nomenon^^. Further away from the hole mass-shell the 
spectral function crosses over to the TL one, Eq. (4). 

We described above the behavior of the spectral func- 
tion close to the mass-shell. The regions of a finite spec- 
tral weight are indicated by the shaded area in Fig. 1. 
Apart from the hole mass-shell, all edges of the spec- 
tral support arc characterized by the power-law behavior 
of the spectral function with momentum-dependent posi- 
tive exponents. This is easy to foresee, if one notices that 
the states responsible for the nonzero spectral weight far 
from the single-particle mass-shell must involve several 
excitations. The main contributions to the values of the 
exponents come from the constraints on the phase space 
available for such excitations. The additional contribu- 
tion to the exponents is due to the interaction between 
the excited particles and holes. This contribution is first- 
order in the interaction potential and has the same ori- 
gin as the the exponents in the X-ray singularity phe- 

23 



the structure of the quasiparticle peak near the particle 
mass-shell. 

In the hole region of the spectrum e < 0, the positive 
dispersion curvature does not smear the TL power-law 
singularity at the mass-shell. The singularity is preserved 
because due to kinematic constrains, the hole mass-shell 
remains to be the edge of the spectral support (see the 
next section). The only effect of the spectrum nonlin- 
earity is the renormalization of the exponent in the hole 
part of the spectral function [cf . Eq. (4)] , 



Akie) oc 



(13) 



where — e <C k"^ /2m and the exponent 7^ is given by 
Eq. (8). Note that the exponent is invariant upon mo- 
mentum inversion k ^ 2pp-\-k and interpolates smoothly 
between the two Fermi points; Eq. (13) is applicable 




FIG. 2 

shell e 



Spectral function in the vicinity of the particle mass- 
si > 0. The sohd hne is a plot of Eq. (7) with 
^^2pp+k = 1.1, Mfc = 0.4, and 7fe given by Eq. (8). For com- 
parison, the dashed line corresponding to the Luttinger liquid 
result Eq. (4) is also shown. 



of the spectral function is developed in Section V. 

The rest of the paper is organized as follows: In Sec- 
tion II we present a qualitative analysis of the problem. 
Perturbative calculation of l/rj, to fourth order in the in- 
teraction strength is carried out in detail in Section III. 
In Section IV we develop a strategy of summing up the 
leading logarithmic corrections while accounting for a fi- 
nite I/t^, and derive Eq. (7). The behavior near the 
edges of the spectral support is discussed in Section V. 
We compare our findings with the results obtained for 
the exactly solvable Calogero-Sutherland model in Sec- 
tion VI. Finally, discussion and outlook are presented in 
Section VII. 



II. QUALITATIVE CONSIDERATIONS 

In this section we discuss the boundaries (kinematic 
edges) of the area in (p, e)-plane where the spectral func- 
tion differs from zero (see Fig. 1). We also provide a 
simple qualitative framework for understanding the be- 
havior of the spectral function near the kinematic edges 
and in the vicinity of the particle mass-shell. 



A. Support of the spectral function 

Using the Lehmann representation, we may express the 
particle contribution to the spectral function as 

^fc(e)«EK/i<^|o)l''^(^-^i/>+^io>)^(^-^i/>+^io>)- 

I/) 

(14) 

Equation (14) may be viewed as the probability of tun- 
neling of a particle with a given momentum pp + k and 
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energy ep + e into a ID system (to be definite, we con- 
sider the right-movers pp + k > 0). The initial state of 
the transition in Eq. (14) is the ground state of the sys- 
tem |0), and the final state is |/), with J]o); P\f) a^nd 
i5|o), E\f^ being the corresponding momenta (relative to 
Pp) and energies (relative to ep), respectively. 

In the absence of interactions, the only possible final 
state is that with a single right-moving particle added 
to the groimd state. Because of the momentum conser- 
vation the added particle must have momentum k (here- 
inafter the single-particle momenta are measured relative 
to the respective Fermi points), |/) = t/jj^'^ |0). Equa- 
tion (14) then yields ex 9{k)6{e - ^j}). With inter- 
actions present, states |/) may contain, in addition, a 
number of particle-hole pairs. This allows Ak{e) to be 
finite away from the single-particle mass-shell e = 
Still, there are regions in (p, e)-plane where Ak{e) = 
due to kinematic constrains. To simplify the analysis of 
these constraints we focus on the low energies, |e| <C ep, 
and small momenta |fc| <C kp. 

The simplest low-energy final state, beyond the single- 
particle one, is the state containing one additional 
particle-hole pair. Given that the total momentum (rel- 
ative to Pp) is small, the additional pair must have a 
small momentum too. Therefore, it has to be located in 
the vicinity of either the right or the left Fermi point. 
These states, see Figs. 3(a) and 3(b), have the form 

l/> = ^kX^l |0> , fci > , fc2 > , fcs < (15) 
and 

I/) = V'f>fcJ< |0) , fci > , fc2 < , fcs > 0, (16) 
respectively. The state (15) has energy 

- £;|o) = v{ki + &2 - fcs) + ^ {kl + kl- kl) . (17) 

Taking into account the momentum conservation k\ + 
k2 — ks = k, we find 

E\f^ - E\o^ =vk + ^[kl + kl- {kl +k2- kf] . (18) 

The constraints on ki, k2, and ks in Eq. (15) guarantee 
that fc > 0. At a given k, the smallest possible value of 
the excitation energy E^f^ — E^q-^ is reached at ki = k2 = 
0, fcs = -k. 

Similar consideration for the state (16) yields (with 
momentum conservation ki + k2 — ks = k taken into 
account) 

^,^-E^o^=-vk+2vki-^ + ^{ki+k2)-^. (19) 

At A; < 0, the lowest energy is reached at fci = ^2 = 0, 
fcs = -k. 

An important observation following from the above 

analysis is that the lowest possible excitation energy cor- 
responds to the final states with all particles at the Fermi 



level, and a single hole with the largest possible absolute 
value of the momentum. It is easy to check that a fi- 
nal state with more than one particle-hole pair excited 
from the ground state still has the lowest energy when 
all the available momentum is "carried" by just a single 
hole. Thus the final states (15) and (16) correspond to 
the kinematic boundary for the particle (e > 0) part 
of the spectral function. Combining Eqs. (18) and (19) 
at kl = k2 = 0, we find 



^fc= v\k\ — k^/2m. 



(20) 



Although we considered here the small-fc domain only, 
Eq. (20) is valid for the entire region —2pp < k < 2pp, 
see Fig 1 . Finding the spectral edges for higher momenta 
|fc| > 2pF, however, requires consideration of final states 
with more than one particle-hole pair. 






FIG. 3: Final states |/), see Eqs. (14) and (21), corresponding 

to the boundaries of the regions whore Ak{e) > 0. The black 
(white) circles stand for extra particles (holes) added to the 
ground state; the dashed line indicates the Fermi level. States 
(a) and (b) correspond to the particle edge (e > 0); states (c) 
and (d) correspond to the hole edge. 

When considering the hole part of the spectrum, we 
again start with the Lehmann representation, 

^fe (e) ^ |(/ 1 Vl" |0) I -P,o) ) ^ -i?|0) ) . 
l/> 

(21) 

Here P|/) is the momentum of the final state relative to 
—pp. In the absence of interactions, the matrix element 
in Eq. (21) is finite only for a single hole excitation with 
momentum fc < 0, and A), oc 0{—k) S{e — ^^). Interaction 
results in a finite Ak{e) at e 7^ Indeed, consider the 
state [see Fig. 3(c)] 

I/) = «C |0) , fci < , fc2 > , fcs < 0. (22) 
In zero order in interaction, the energy of this state is 

E\f^-E\o)=vi-ki+k2-k3)-^ {kl + kl-kl). (23) 

To find the spectral edge [i.e., the largest possible value 

of e in Eq. (21) at a given fc], we look for the lowest pos- 
sible energy of the final state E^f^ in Eq. (23). Given 
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the momentum conservation, ki + k2 — = k, this hmit 
is reached at ki = k, k2 = k^ = for A: < 0, and at 
k2 = k, ki = k^ = for A: > 0. In other words, the 
lowest-energy final state (hence the highest possible e) 
coincides with the energy of a single hole, see Figs. 3(c) 
and 3(d). Consideration of final states with more than 
one particle-hole pair excited does not change this conclu- 
sion. Accordingly, the kinematic boundary for the hole 
part of the spectral function at small k is given by 

0^ = -v\k\ + k'^/2m. (24) 

Just as it is the case for e > 0, Eq. (24) is valid for 
I A; I < 2pf; finding the kinematic boundary for higher 
momenta requires consideration of final states with more 
excitations. 

The support of the spectral function in (p, e)-plane is 
illustrated in Fig. 1. 



B. Spectral function nesir the edges of the support 

Near the edges of support in {p, e)-planc. the spectral 
function displays a power-law dependence on the distance 
to the edge. The behavior of Afc(e) near the hole mass- 
shell, 

^fc(e)<x(a-e)^^'"'^(a-e), -2pF<k<0, (25) 

is very similar to that in the TL model. The deviation 
of the exponent in Eq. (25) from (—1) is quadratic in the 
interaction and positive, just like it is in the TL model. 

Adding an extra hole to an interacting system is ac- 
companied by the creation of "soft" particle-hole pairs; 
we considered an example of such three-particle fi- 
nal state in the previous subsection, see Eq. (22) and 
Fig. 3(c). Excitation of multiple low-energy particle-hole 
pairs leads to the orthogonality catastrophe^^, thus the 
correction to (—1) in the exponent is positive and is pro- 
portional to 7^. The main difference compared with the 
TL model is that now the extra hole interacts with the 
soft pairs residing near both Fermi points, while in the 
TL model only pairs on the opposite branch of the spec- 
trum contribute to the singularity. 

When the edge of the spectral support no longer co- 
incides with the mass-shell, which is always the case for 
particles (e > 0) and at < A: < 2pp also for holes 
(e < 0), the behavior of the spectral function near the 
edge is different from that described by Eq. (25). We 
find that the corresponding exponents are positive and 
finite even in the limit 7^ ^ 0. 

To be definite, let us consider the particle edge e 

at A: < 0. The spectral function differs from zero in 
the vicinity of the edge due to three-particle final states 
I/) = V'f>feJ< see Eq. (16) and Fig. 3(b). If we 
replace the matrix element in Eq. (14) by a constant, 
then the dependence of (e) on e — can be found by 



power counting, 

Afe(e) cx J dkidk2dk^ 5{ki + A'2 — k^ — k) 

(X {e-^e{e~^). (26) 

Eq. (26) corresponds to the final state |/) with two 
particles close to the Fermi points ±.pF and a "deep" 
hole with momentum (relative to —pp) approaching k, 
see Fig. 3(b). Interaction of the two particles near the 
Fermi level with the hole and with each other leads to 
a logarithmic renormalization of the matrix element in 
Eq. (14). Similar to the X-ray edge singularity prob- 
lem^^, the renormalization occurs already in the first or- 
der in the interaction strength and leads to 

Ak{e) oc (e - ^^'^-i^>^-i^2pp+k+'2iJ.2pp ^ _2pp < fc < 0. 

(27) 

For a short range interaction and small momenta |A;| <C 
Pf we have ^,k <C At2pj? = 270, and Eq. (27) simplifies to 

^fe(e)cx(e-a)'+'^°. 

For a linear spectrum (TL model), the left-moving parti- 
cle and hole have the same velocity, destroying the core- 
hole effect leading to Eq. (27). 




(28) (29) (7) 



A:<0 . 








^f- 





(25) (27) 



FIG. 4: Dependence of the spectral function on e at a fixed 

k — p — Pf with A; > (upper panel) and A- < (lower panel) . 

The numbers in brackets refer to the equations describing the 

corresponding asymptotes of Ak{e). At large |e— ^fc|, |e— ^fc| ^ 

k'^/m (this region is not shown in the figure) the conventional 

2 ^ — 1 

Luttinger liquid behavior ^fe(e) oc \e + vk\''° \e — vkl''" is 
restored. 

The consideration of spectral function near the hole 



edge e— >^featfc>Ois very similar, yielding 

Ak{e) (X (i;^- e)i-2*'2pF-;*fc+/i2j,p-fc ^ Q < fc < 2pF. 

(28) 

The signs of two terms in the exponent here are different 
from those in Eq. (27) because of the difference in the 
structure of the relevant final states. In the three-body 
sector, for example, the final state consists of one particle 
and two holes (rather than two particles and one hole as 
in the case of Eq. (27) above), see Fig. 3(c). For a short- 
range interaction and small momenta, Eq. (28) simplifies 
to 

^fe(e)oc(a-e)'-'^°- 

Finally, we discuss the particle edge e ^ at /c > 0. 
A new element here is that the relevant soft excitations 
reside on the same branch of the spectrum. When e is 
close to these excitations tend to occupy almost iden- 
tical right-moving states with k 0. In the specific 
example of the final state shown in Fig. 3(a). the mo- 
menta of two particles difi^er by (e — ^k)/v. This re- 
sults in a suppression of the matrix element in Eq. (14), 
I (/I V'^^|0) \ oc e — ■ A proper modification of Eq. (26) 
then leads to Afe(g) cx (e — £^k)^. Accounting for the 
particle-hole interaction in the final state yields 



^fc(e)cx(e-a)'~'^N 0<k<2pF. 



(29) 



At k = 2pF (i.e. p = 3pf) the spectral edge ^i- touches 
zero, £_2pF = 0. Close to this point the exponent can 
be evaluated using the conventional bosonization tech- 
nique^^. It yields 3 — 470 for the exponent in agreement 
with Eq. (29) (indeed, iJ,2pF = 27o)- 

The dependence of the spectral function on e at a fixed 
k is sketched in Fig. 4. 



C. Spectral function near the particle mass-shell 

For a nonlinear spectrum (6) the particle mass-shell 
e = > lies above the lower edge of the spectral sup- 
port, ^k= ^k' ~ k'^/fn. Far away from the mass-shell the 
difference between and is not important, and at 
e — » k'^/m the spectral function Ak{e) approaches 
the TL form Eq. (4). In the vicinity of the mass-shell, 
however, ^fe(e) imdergoes a dramatic change. Indeed, 
since £ = £,k no"^ li^s within a broader continuum, the 
quasiparticle relaxation (decay) is no longer prohibited 
by the conservation laws. As a result, in a parametri- 
cally wide region (indicated by dark gray in Fig. 1) the 
quasiparticle peak in Ak{e) acquires a Fermi-liquid-like 
Lorentzian shape, see Eq. (11). 

It is instructive to discuss the origin of a finite relax- 
ation rate in a ID Fermi system. First of all, unlike in 
higher dimensions, relaxation in ID can not occur via 
two-particle collisions. Indeed, conservation of the mo- 
mentum and energy allows at most a permutation of the 
momenta of two colliding particles. 




-Vf+K 




PF+k' 



-PF+k 



FIG. 5: The initial (left panel) and the final (right panel) 
states of a three-particle scattering process that leads to a 
quasiparticle relaxation in one dimension. 



Scattering processes that result in a redistribution of 
the momenta and thus potentially lead to a finite relax- 
ation rate must involve at least three particles. In such 
three-bodycollision three particles with momenta k, ku, 
and kL (relative to ±pf) in the initial state \i) end up 
in a final state |/) with different momenta k' , k'j^, and 
k'^, see Fig. 5. For a generic interaction^^ the transition 
\i) I/) has a nonvanishing amplitude A. 

In order to estimate the relaxation rate of an extra 
right-moving particle; with momentum k due to three- 
body collisions, we note that the single-particle states 
kn, kL in the initial state of the transition \i) are below 
the Fermi level, while all three single-particle states in 
the final state |/) are above it. Applying now the Golden 
Rule, we find 

— oc / dk'dk'ndkLf dk'^dknlAl^ (30) 

Tk Jo J-00 

X S[{k + kR + kL) - {k' + k'ji + ki)] 



xS 



where A is the three-body collision amplitude introduced 
above, and the (5-functions express the energy and mo- 
mentum conservation. 

In writing Eq. (30) we took into account that for 
k <^ Pf the conservation laws cannot be satisfied un- 
less the collision involves both the right- and the left- 
moving particles. Further analysis shows that the con- 
servation laws allow a small (< k"^/ m) momentum trans- 
fer to the left-movers. Such solution can be found by 
iterations. To zero order in kj^ — k'j^, the momentum 
conservation gives k — k' = k'j^ — kj^. The energy re- 
leased in the collision of two right-moving particles then is 

+ ^fc" ~ik' ~^k' ^ ^'^/ This energy is transferred to 
the left-movers, fh — f f ^ k'^/m, which corresponds to 

the momentum transfer kj^ — k'j^ < k'^/mv <C k. Accord- 
ingly, the energy and momentum conservation restrict 
the range of the momenta contributing to the integral in 
Eq. (30) to 



k',k'j,,\kj,\<k, kL,\k'L\<k^/mv. 



(31) 



The (5-functions in Eq. (30) remove the integrations 
over k'^ and k'^. With the phase space constraints (31) 
taken into account, the remaining integrations then yield 

the estimate 



2lA 



l/Tk oc \AYk 



(32) 
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For a weak gcncric^^ interaction, the nonvanishing 3- 
particle collision amplitude A appears in the second order 
in the interaction strength. More careful consideration 
(see Section III) which accounts for the indistinguisha- 
bility of the two right-moving particles participating in 
the collision results in the estimate 



AocVo{Vo-Vk)ock'', 



(33) 



cf. Eq. (10). Accordingly, the quasiparticle decay rate 
scales with k as 1/rfc oc A:*. 



clarity, we consider the decay of a right-moving particle 
(^^ > 0). In the presence of the spectrum nonlinearity 
ImEfc acquires an additional contribution in the second 
order in the intrabranch interaction V. This contribu- 
tion comes from the scattering processes with two parti- 
cles and one hole in the final state, all three on the right 
moving branch. The Golden-rule expression for ImEfe 
reads 



■I-5:f(e) = |E 

fci ,k2 



S{e-C^,-^g+^g), (36) 



III. PERTURBATION THEORY 

In this section we evaluate the imaginary part of 
the self-energy ImEp(e) perturbatively in the interaction 
strength. Specifically, we focus on the near vicinity of 
the mass-shell e w ^p. The Fermi liquid theory predicts 
a nonzero self-energy already in the second order in the 
interaction strength, see Eq. (2). We show below that in 
one dimension the second-order contribution to ImEp(e) 
vanishes on the mass-shell, even if the curvature of the 
dispersion relation Eq. (6) is taken into account. A finite 
quasiparticle decay rate appears only in the fourth order, 
and only if ^ 0. 

We describe interacting spinless fermions by the Hamil- 
tonian 

a,k q^O 

(34) 

where a = R,L, the dispersion relation is given by 

Eq. (6), and = X^fcV'fclqV'fc is the Fourier component 
of the density operator. We found it convenient to distin- 
guish between the interbranch and the intrabranch inter- 
action potentials (denoted by Uq and Vq, respectively)^^. 
However, for the sake of brevity, we will set Uq = Vq in 
the results. 



A. Second order 

Evaluation of the self-energy in the second order of per- 
turbation theory at |e| > |^p| parallels the corresponding 
calculation in the problem with linear dispersion relation 
(Tomonaga-Luttinger model). The only finite contribu- 
tion comes from the interbranch interaction U, and the 
result reads 



where k — p—pp and ki + k2 — k^ = k due to momentum 
conservation. The Fermi statistics dictates that ki,k2 > 
and ^3 < . The amplitude of the scattering process is 
given by 



(37) 



The conservation laws arc satisfied only in the limited 
energy range, ^ < e < S,^ with ^ = - P/m, see 
Eq. (20). We defer the discussion of the spectral function 
near the lower edge of the spectrum e = to Section V, 
and focus here on the immediate vicinity of the mass- 
shell, |e - < k'^/2m. In this Hmit Eq. (36) yields 



■ImSf 



(e) 



l^k (fR 
-J \Sk 



(38) 



with /ifc defined in Eq. (9) and given by Eq. (6). The 
real part of the self-energy is given by 



-ReEf (e) 



Ml 
47r 



k'^/m 



(39) 



According to Eq. (38), the intrabranch interaction 
along with the positive curvature of the dispersion re- 
lation (1/m > 0) results in a finite ImE^^^ at e < 
Eq. (38) complements Eq. (35), familiar from the con- 

ventional TL theory. Note that in both cases Im EJ^ ' (e) 
vanishes on the mass-shell due to the phase space con- 
straints. 

For the particle branch of spectrum (e > 0), the 

kinematic restrictions on the phase space are lifted in 
higher-order processes. The simplest process leading to 
a finite self-energy on the mass-shell is the second-order 
process that results in the creation of two particle-hole 
pairs in the final state. 



1 



In- 



■*^(kl-l?pl) 



^^^^ 

with 7o = vUq/2. Vanishing of the imaginary part of the 
self-energy below the mass-shell is due to kinematic con- 
straints: the phase space available for scattering process 
vanishes in the limit e — » ^p. 

For |e| < l^pl real decay processes are allowed by con- 
servation laws only if the dispersion is nonlinear. For 



B. Fourth order 

There are two kinds of the fourth order processes. The 
first one leads to a logarithmic correction to the scatter- 
ing amplitude via virtual creation of a particle-hole pair. 
Similar logarithmic corrections appear in higher orders of 
perturbation theory as well. Summation of the leading 
logarithmic contributions in all orders (the corresponding 
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procedure is described in Section IV) results in a power- 
law behavior of the spectral function. This is very simi- 
lar to the TL model, although with the exponent slightly 
modified due to the dispersion nonlinearity. 

A different kind of second-order processes, leading to 
the finite on-shell value of ImSfe, involve creation of five 



quasiparticles in the final state: a particle on the right- 
moving branch, and two particle-hole pairs. Kinematic 
considerations show that the two pairs must be excited 
on the opposite branches of the spectrum in order to yield 
ImE;j 7^ on the shell. The corresponding contribution 
to the self-energy is then given by 



-lmEW(e) 



TT 



2 £ 

g,Q,fei ,fe2 



iki ik^-Q 



(40) 



r 



where the summation range is limited by the Pauli prin- 
ciple constraints 

k - q + Q, ki, k2 > ki - q, k2 - Q < 0, 
and the amplitude A^^^ consists of two contributions, 

A^ ^ = •^uir + -^uv • i'^^) 



These two contributions correspond to two possible ways 
the desired final state can be reached in a second-order 
process: The first one, A^jj oc U'^, arises solely due to 
the interbranch interaction U, see Fig. 6. The second 
one, A^y oc UV , involves both the interbranch {U) and 
the intrabranch {V) interactions, see Fig. 7. The corre- 
sponding analytical expressions read 



A 



(2) 
UV 



+ 



cR 

^ki—q 



^ki + ^k2 ^k2-q 



Uk-ki+gUk-ki 



Uk-ki+qUk-ki 



^k-q+Q Sfci_gT^?fc2-Q 

UQp^q-Q - Vk-ki+q) 
fR cR cL I cL 

?fel ?/ci-Q ?fc2 ?/S2-Q 

UQ{Vq-Vk-k,) 



cL cR cR 1 cL 

^k2+k-ki '^ki-q '^k-q+Q ^ ^k2 



(42) 



2—Q+ki—k 



+ 



UQ{Vq.Q-Vk-k,) 



'^ki-q 



cR -L (Ti 



?fe2-Q 



+ 



fR 



q+Q 



^k-q 



^k2 + ^k2 



+ 



pR 



Vq{V, - Vk-k,+Q) 

— i:R j-CR — fR 
q+Q ^k-L-q^^k-L '^k+Q 



(43) 



(A more general but rather cumbersome expression for 
the amplitude free from the simplifying assumption [21] 
is derived in Ref. [25].) 

One immediately notices that A\jy vanishes for the 
momentum-independent intrabranch interaction Vq = 

Vq. Although it is not obvious, the amplitude Auy van- 
ishes'^ for momentum-independent interbranch interac- 
tions Uq = Uo and for a strictly parabolic dispersion re- 
lation. 

To proceed further, we assume that both inter- and 

intrabranch interactions are symmetric analytic functions 
of the transferred momentum; at small momenta 



Uq^Uo + ^ U'^q^ , 



(44) 



and at large q the potentials vanish suSiciently fast'^. 
We also neglect cubic and higher-order terms in the dis- 
persion relation Eq. (6). 



We will evaluate Eq. (40) for e = In this limit the 

typical momentum of the right-moving particles in the 
final states contributing to ImS^'*^ is of the order of k 
and the velocity variation is k/m. The energy gain 
due to the production of the right-moving particle-hole 
pairs is therefore ~ k{k/m). Conservation of energy then 
yields the estimate for the typical momenta of the left- 
moving particle and hole in the final state, 

k^ <^ Q <^ k^/mv <C k. 

Carrying out the summation over ki in Eq. (40), we ob- 
tain 

q,Q 0<k2<Q 

where the summations over q and Q are restricted to the 
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J? 



U 



u< 



R 



K-Q 
k-g+Q 




(a) 



(b) 



FIG. 6: Second order in U contributions to the amplitude 



A 



(2) 



and (b) correspond, respectively, to the first and 



second terms in Eq. (42) . Two more contributions (not shown 
in the figure) correspond to the replacement fei <-> fe — g + Q 
in the final states; these are the third and the fourth terms in 
Eq. (42). 





k-q+Q 




K-1 

R k-q^Q 



FIG. 7: Second order contributions to the amplitude 
(a)-(d) correspond to the first terms in the numerators of 
the four contributions in Eq. (43). Four more contributions 
(not shown in the figure) correspond to the replacement fci <-» 
A; - g + Q in (a)-(d). 



domain 



< < fc, 



{k-q)q 
2mv 



<Q< 



kq 

2mv' 



(46) 



The amplitude A^'^\ see Eqs. (41)-(43), simplifies con- 
siderably in the limit k <C mv. Keeping only linear in 



Uq and Vq contributions, we find 



- {2mvQy 
Smv'^q'^ 



(47) 



The high powers of the momenta here as well as the factor 
1/m resulted from delicate cancellations among various 
contributions to the amplitude. Since the amplitude (47) 
is independent of fc2, the integration over this variable in 
Eq. (40) brings about a factor Q/2tt. The remaining 



integration over q and Q is restricted to the domain D 
defined in Eq. (46), 

Jd 477^ Snq 
The integration here is straightforward and yields 

_^^,,,^„,^MiHa±i£p!. (48) 



5 (327rmt;2)3 



The fac- 



Note that ImE^^(^^) scales with k as i 
tor /c^ here originates in the amplitude: according to 
Eq. (47), its typical value is A^"^^ oc P. The remain- 
ing factor comes from the integration over Q (re- 
call that Q oc k^). Using now /c^Vq" « 2{Vk — Vq), see 
Eq. (44), and setting Uk = Vfe, we arrive at Eq. (10) with 

l/2Tfe = — ImS^'*^ (^j^) [here we took into account that 

ImSf (a)=0]. 

The above derivation can be extended to e 7^ The 
self-energy varies with e on the scale of the order of 
k'^/2m, and it vanishes identically for e < ^fc, as expected 
from the kinematic considerations of Section II. 

According to the developed perturbation theory the 
spectral function in the vicinity of the particle mass-shell 
takes the form 



1 



ImS,(c) 



^ (e-4«)^ + [ImS,(e) 



(49) 



where Sfc(e) = ^i^'(f) + ^''k \^-) the corresponding 
contributions arc given by Eqs. (35), (38), and (48). As 
a result, the spectral function in the energy interval 



around the mass-shell is a Lorontzian with the width 



(2) 

l/2rfe. Outside this interval ImS^^(e 
and the spectral function is described by a power law 
Afe(e) ^ |e — Cfe^l"^, with the exponent equal to (--1). 
The perturbation theory developed so far does not take 
into account the logarithmic renormalization of the scat- 
tering amplitudes. In Section IV we show that summa- 
tion of the leading logarithmic corrections leads to the 
interaction-dependent correction to the exponent. The 
height of the Lorentzian peak Eq. (49), is also renormal- 
ized, see Eq. (11) above. 



> ImE 



(4)/ 



IV. LEADING LOGARITHMIC CORRECTIONS 

It is easy to see that excitation of virtual particle-hole 
pairs leads to the logarithmically divergent contributions 
in perturbation theory, in a very much the same way as in 
the TL model. The real part of the self-energy acquires 
such logarithmic corrections already in the second order 
in interaction, see Eq. (35). The logarithmic terms in 
the imaginary part of the self-energy appear in the fourth 
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order. [These contributions vanish on the mass-shell, and 
thus do not affect the vahdity of Eq. (48)] . 

In this section we develop a procedure for the summa- 
tion of the leading logarithmic corrections to the Green 
function Gfc(e) in the presence of nonlinear terms in the 
dispersion relation Eq. (6). The nonlinearity is not im- 
portant as long as |e — ^ k'^/2m. Close to the mass- 
shell, however, the behavior of the spectral function de- 
viates significantly from that in the conventional Lut- 
tinger liquid. For a positive curvature [i.e. for 1/m > 
in Eq. (6)], deviations are the strongest for the particle 
(e > 0) excitations with fc > 0; we have already seen that 
the on-shell excitations acquire a finite lifetime Tk- On 
the other hand, the behavior of the hole branch (e < 0) is 
qualitatively similar to that of a Luttinger liquid. We will 
concentrate here on the properties of the particle branch 
at |e — Cfcj ^ k'^/m, deferring the discussion of the hole 
region of the spectrum till the end of this Section. 



A. Vicinity of the particle mass-shell: e 



fc > 



In order to account simultaneously for both the loga- 
rithmic renormalization and the finite quasiparticle life- 
time, we notice that the relevant energy scales form a 
well-defined hierarchy 



(50) 



The logarithmic corrections originate from almost the en- 
tire energy band, i.e. from the states with energies in the 
range (|e — ^^|, mt;^). On the other hand, the finite life- 
time 1 /Tk originates in the decay of a particle into particle 
and hole states within much narrower strip of energies of 
the width of the order of 



(a) 



Ak 




Pf 
Xk 



V 



d 



d 



FIG. 8: (a) Subbands in the momentum space remaining after 
performing the RG transformation described in step (i) in the 
text, (b) Subbands of the effective model described by the 
Hamiltonian (58). The state with fc > under consideration 
belongs to d-subband. 

Our strategy for evaluation of the Green function will 
be as follows: 

(i) First we will take into account virtual transitions 
to states with relatively high energies, corresponding to 

the momenta p in the range Afc ^ \p — Pf\ ^ Pf for 
the right- movers and Ak^/mv < + ^ Pf for the 



left-movers (here A 3> 1). This step is similar to the 
standard renormalization group (RG) calculation for the 
Tomonaga-Luttinger modeP''. As a result of "integrat- 
ing out" the high-energy states, Gfe(e) acquires a multi- 
plicative factor. In addition, some irrelevant terms are 
generated in the Hamiltonian; these terms do not affect 
the logarithmic renormalization of Gk{e), but contribute 
to the decay rate l/rfc which is evaluated on the lat- 
ter stage. The multiplicative RG terminates when the 
energy bands have been reduced to strips of the right- 
moving states with \p — Pf\ ^ Afc, and the left-movers 
with Ip + pfI Ak'^/mv, see Fig. 8(a). 

(ii) In the next step, we separate the remaining states 
in three groups. The first two (subbands r and I) cor- 
respond to two segments, \p — Pf\ ^ Afc and |p ^ 
Xk^/mv on the right- and left-moving branches, respec- 
tively; hereinafter A ^ 1. The rest of the states form 
subbands d, see Fig. 8(b). The state of interest fc belongs 
to the right-moving branch of rf-subband. This state ac- 
quires a finite lifetime Tk, which manifests itself in the 
shifting of the pole in Gfe(e) by i/2Tk off the real axis. 
Importantly, the dominant contribution to the finite life- 
time comes from the decay of the state of interest into 
other states within d-band. Although the subbands r and 
I make a negligible contribution to the decay rate l/r^, 
the density fluctuations in these subbands induce a slowly 
varying [on the timescale of the order of (fc^/2m)~^] fluc- 
tuating field that affects the dynamics of the high-energy 
(e ~ d-particle. Our formalism accounts for both 
these fluctuations and for the finite decay rate l/2rfc, 
and leads to Eq. (7). 

Now we sketch the implementation of the two steps 
described above. (Some technical details are relegated to 
Appendixes A and B.) 

Step (i) follows, with a small modification, the con- 
ventional RG procedure employed in the theory of 
Tomonaga-Luttinger modeP^. We introduce a reduced 
space of single- particle states with \p ~ pf\ < fcfi for 
the right-movers and |p + pf| < for the left- movers 
and consider the Green function, G;j(e; fcz,, fcfl), defined 
in the reduced band. Following the idea of the multiplica- 
tive RG^^, we consider the transformation of the Green 
function associated with the reduction of the bandwidth 
(fci — > fc^, fc/j ^ fc^), and cast it in the form 



Gk{e;klk'j,)=z( ) Gk{e;kL,kn) 



(51) 



Next, we use the results of the Section III A to evalu- 
ate the Green function (51) perturbatively in the second 
order in the interaction potential. This leads to an ap- 
proximate expression for z, 

where Afc < fc^ < pf- The last term on the right- 
hand side here expresses the fact that the renormaliza- 
tion of Gfe(e) for a right-mover results from its inter- 
action with the left-movers; this is why the function z 
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depends only on a single argument fc^/fc^. Eq. (52) is 
valid up to the seeond order in 70, and also assumes 
that 7oln(fc^/fci) ^ 1. In order to find z{k'^/kL), we 
supplement Eq. (52) by the requirement that the scaling 
function has a multiplicative property^''. That leads to 

(53) 

In order to keep the state of interst k under consid- 
eration it must lie within the reduced band of right- 
movers. Accordingly, the scaling of the right band must 
stop at A;^ « Afc. However, the bandwidth of the left- 
movers can be reduced even further, all the way down to 
k'j^ K, hk^jmv. The latter scale represents the width of 
the band to which the states involved in real transitions 
associated with the formation of 1/rfc belong, see Sec- 
tion IIIB. (Note that according to Eq. (35), the scaling 
exponent in the interval Kk'^/mv < k'j^ < Ak is 7q rather 
than 27q). 

Using now Eqs. (51) and (53), we find 

fAkvY-^'^fAe/my'" ( AP ^ \ 

(54) 

(we used mv ~ here). Note that the nonlinearity of 
the single-particle spectrum does not affect the form of 
Eqs. (52) and (54) as long as A > 1. 

However, because of the nonlinearity, reduction of the 
bandwidth {kp — * Ak'^/pp) generates an additional in- 
teraction term 

= ^M^f>fI,_Q<t_^Vfe.V'^<-, . (55) 

The matrix element M here accounts for the contribu- 
tions to the total transition amplitude A^^^ of the virtual 
states outside the reduced band of the left- movers. These 
states contribute to Ajjjj , but not to Ajjy . Therefore one 

(2) 

finds M fa {l/2)Aijlj (here 1/2 is a combinatorial factor). 
Since all the momenta in Eq. (55) belong to the reduced 
band, we have 

as it is clear from the analysis leading to Eq. (47). Us- 
ing q ^ k and Q ~ k'^/pp, we estimate M ~ fc^/mw^, 
i.e., M vanishes for a linear spectrum. The generated 
interaction Eq. (55) is to be added to the Hamiltonian 
Eq. (34), where all the operators now act within the band 
of reduced width. Along with the interaction terms al- 
ready present in Eq. (34) , V\ contributes to the inelastic 
scattering processes giving rise to a finite quasiparticle 
relaxation rate 1/rfe. 

We now proceed to step (ii) of the program outlined 
above, i.e., evaluation of the Green function in the re- 
duced band. To this end, we separate the states in the 
reduced band in two narrow subbands: the i-band with 



1^*1 ^ Xk'^/pF around the left Fermi point, and the r- 
band with < Afc around the right Fermi point. The 
remaining states form subband d, see Fig. 8(b). Param- 
eter A here satisfies the conditions 

7^1n:^<i;l, A«:l. (57) 
A 

(Note that at 70 <C 1 the two conditions (57) can be 
satisfied simultaneously.) The states within the subbands 
r and I produce additional logarithmic renormalization of 
the Green function. The first condition in Eq. (57) allows 
us to disregard such corrections originating in d-subband. 
The second condition in Eq. (57) ensures that accounting 
for the inelastic scattering within d-subband reproduces 
the result of Section HI for the relaxation rate l/r^. 

We write the Hamiltonian H+Va acting in the reduced 
subspace r,l,d (see Fig 8) as 

H + VA = Ha + Hrl+Hd-rl. (58) 

The first and second terms in the right-hand side here 
account for the states in "high-energy" subband d and 
the low-energy subbands r and I, 

Hd =Pd{H + VA)Pd, Hrl = Prl {H + VyPrl , 

where Pd and Pri are projectors onto the corresponding 
states. The third term in Eq. (58) describes the interac- 
tions, 

Hd-rl = Pd{H + VA)Prl + Prl{H + Va)^^. 

Using Hamiltonian Eq. (58), we evaluate the Green func- 
tion of a particle in the subband d. 

We start the evaluation of Gh (e;^^,Ak) with ac- 

counting for the interaction terms in Hd- Evaluation of 
the self-energy follows the perturbative analysis of Sec- 
tion HI. The only difference is that the second-order 
contribution to ImS^, see Eq. (38), is absent provided 

that |e — ^1^1 < Xk'^/mv. Indeed, a finite ImSj,^^ re- 
quires excitation of a particle-hole pair in the Z-subband. 
Such processes are absent in Hd but are included in 
Hd-ri discussed below. In addition, the second-order 
contribution to the real part of the self-energy. Re ^ 
7Q(e — ^1*^) ln(A/A), can be also neglected as it is small 
compared with [e — ^^[. Unlike in the second order, the 
low-energy states {I and r subbands) do not play any spe- 
cial role in the fourth order calculation, see Section HI. 
Hence, Hd alone is sufficient in order to reproduce the re- 
laxation rate 1 / 2Tk (note that the term Va is important 
for this calculation). Thus, neglecting Hd-ri, we find for 
the retarded Green function 

Gt(e;—,Ak]= . ,^ . (59) 

mv' J e-^^ + i/2Tk ^ ' 

Next, we account for the effects of Hd-ri interaction 
on the Green function of a particle in the d-subband. For 
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this calculation V\ is irrelevant. Moreover, the interac- 
tion between r and I subbands does not contribute to 
the Green function in the leading logarithmic approxi- 
mation (this interaction leads merely to the higher-order 
corrections to the exponents). We thus neglect both these 
terms and write Hri and Hd-ri as 

Hrl = E E4"V'fV',« (60) 
Oi=r,l k 

Hd-ri = (Vo-Vk)^ E py-, (61) 

\q\<Xk 

+ {Uo-U2p,+k)2 E 

\q\<Xk^/mv 

I 



where the form of the interaction potentials Vq — Vk and 
Ua — U2pp+k originates in the reduction of the general 
interaction to the density-density form. Because of the 
strict limitation on the wavelengths of the density fluc- 
tuations, the subbands r and I do not contribute to I/t^. 
Therefore the spectrum within these subbands can be 
linearized. Since both r and I subbands contain the 
Fermi level, the excitations within these subbands can 
be described using the conventional bosonization tech- 
nique. Upon introducing the bosonic fields </5^'' such that 
Pg'' = iqifig'^ /27r, we represent Hri and H^-h in the form 



^- = 7^E vi^^r+T^ E -^^k^r (62) 



|q|<AA; \q\<Xk'^/ rnv 



Hd-ri = E 1P>-<i-^ iTzf E IPW-,- (63) 



27rL ^ 2ttL 

|g|<Afe \q\<\k'^/mv 



In the second order of perturbation theory the interac- 
tions in Hd can be neglected. The self-energy of rf-particle 
is then given by the diagram shown in Fig. 9(a). Its eval- 
uation reproduces the logarithmic divergence at e — > 
see Eq. (39). 




FIG. 9: (a) lowest-order logarithmic correction to the d- 
particle propagator (solid line); the dashed line stands for 
r{l) bosons, (b) example of renormalization of the d-particlc- 
boson vortex; wavy lines represent interaction within the d- 
subband. (c) one of the contributions to the decay rate of a 
d-particle. 

To go beyond the lowest order, we need to take into ac- 
count simultaneously both the interactions present in Hd 
as well as the interactions of d-particles with the bosons. 
Since we are interested in the limit e — — > 0, it is suf- 
ficient to consider only the most divergent in this limit 
contributions. 

Clearly, accounting for both types of interactions pro- 
duces two kinds of contributions. First, interactions 
within the d-subband renormalize the particle-boson in- 



teraction vertex. Vertex corrections siich as that in 
Fig. 9(b) can be shown to be small (see Appendix A), 
and we neglect them hereinafter. Second, the interactions 
"dress" the bare d-particle Green function, see Fig. 9(c), 
replacing it by that given in Eq. (59). Obviously, the ef- 
fect of the dressing is to cut the logarithmic divergencies 
off at |e-Cf| - 1/Tfe. 

The task of evaluating contributions such as that 
shown in Fig. 9(c) is greatly simplified by the fact that 
particle-boson interaction is associated with rather small 
(< \k) transferred momentum. This allows us to lin- 
earize the rf-particle spectrum about 

Cl' = ef+fe' ~ + vdk', Vd = v + k/m. 

It is convenient to write the d-particle Green function 
in the real space-time representation, 

Gk{e;^ ,h.k^ = J dx J dte''*Gi{x), x= {x,t). 

In the absence of particle-boson interaction Gf. satisfies 
the equation 

01 +-4-^^'+^) ^^(-)=^(-)- 

With vertex corrections neglected, the effect of the 
bosonic fields is merely to induce a slowly varying in space 
and time potential in which the d-particle moves. This 
fluctuating potential can be written as 

nx, t) = — d^ifi"^ + d^ip' . (64) 
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The retarded Green function of (i-particle in the presence 
of the fluctuating potential G^(x|(/)) satisfies the equation 



dt 



dx 



(65) 

and Gfc(x) is obtained by averaging G^(x|^) over the 
gaussian fluctuations of the field 4>, 

Here {■ ■ ■)4> denotes the time-ordered averaging over the 
slowly varying field (j), see Eq. (64). Note that Gf in 
Eq. (65) describes propagation of a single rf-particlc in 
an empty band, hence the corresponding retarded Green 
function coincides with a time ordered one. 
Carrying out the calculations, we find 

Gi (x|</,) = -ie{t) S{x - Vdt) e-'^^*-'/^^' e'I^W-^(°)l 

(66) 



with 



^iqx—iujt . 



uJ-Vdq 



(67) 



the summation over q here is restricted to < AA: for 
r-d interactions and to |g| < Xk'^/mv for l-d interactions. 
Note that because 9 enters Eq. (66) only in the combi- 
nation 0{vdt, t) — 0(0, 0), the pole at w = v^q in Eq. (67) 
does not show up in Eq. (66). 



With the dynamics of the bosonic fields described by 
a quadratic Hamiltonian Eq. (62), the averaging in Eqs. 
(66) and (67) is straightforward. It yields 



G^(x) = -ie{t) 5{x - Vdt) e 



-if(x) 



with 
if(x) 



^(Wx)-^(o))-^ 



(68) 



(69) 



/ 



duj ((?;)(q)(/)(-q))^ 



27r {u>- Vd qY 



Here 



(0(q)0(-q))^ 

= i{UQ - U2p^+kf tf(q) + i{Vo - Vkf n'-(q),(70) 

and 

Il^'\q) = -tip'-'^p^'')^ = iq/2w){±L0 - vq + tOsgnqy^ 

arc density-density correlation functions for r and I 
bosons. Substitution of Eq. (70) into Eq. (69) and then 
into Eq. (68) yields 



G^(x) = -ie{t)S{x - Vdt) e-'«'=*-*/2^'= exp J -^^^E^\n 



1 + I ix + vt) 

mv 



- ^ In [1 - iXk{x - vt)] 



(71) 



The Green function in the reduced band is given by the Fourier transform of this result (see Appendix B for the 
details of the calculation). 



Xk'^/m 



+ ■ 



AjI \ Xk'^/m 



Ik 



(72) 



Using now Eq. (54) and taking into account conditions Eq. (57), we finally arrive at the expression for the spectral 
function [the abbreviated versions of it were given in Eq. (7) and Eq. (11)] 



1 /f7?\7o / i,2\ 7o-7fc f ,,2 



-1 



I 



Tk \€-i^ + i/2n 



(73) 



The validity of the result (73) is restricted to the 
vicinity of the particle mass-shell, |e — ^ Xk'^/2m. 
It is clear however that well above the mass-shell, at 
e — ^ k^/2m, the dispersion nonlincarity has no ef- 
fect and the spectral function crosses over to the con- 
ventional TL expression Eq. (4). On the other hand. 



e ^ £,k = — k'^ /m, below the mass-shell represents the 
kinematic edge of the spectrum, see Section II. Evalua- 
tion of the spectral function in the immediate vicinity of 
the edge requires a special consideration, and is discussed 
in Section V below. 
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B. Vicinity of the hole mass- shell: k<0 

The behavior of the spectral function near the hole 
mass-shell (e < 0) is nmch simpler because the 

decay of a hole is prohibited by the energy and momen- 
tum conservation laws. The spectral function can be 
evaluated by following closely the route outlined in Sec- 
tion IV A apart from two important modifications: First, 
one needs to replace I/t^ ^ -l-zO in Eq. (65) and its so- 
lution. Second, unlike d-particle, d-hole has a velocity Vd 
which is smaller than v; the difference is due to the posi- 
tive curvature of the dispersion relation. This affects the 
analytical properties of Eq. (71). Indeed, at a; = Vdt the 
factors X + vt and x — vt have the same sign for Vd < v. 
This in turn results in Af^{e) = at e > , which agrees 
with the kinematic constraints on the hole part of the 
spectral function, see Section II. Below the mass-shell, 
at < - e <C ky2m, we find 



Gkie) 



k'^/m J 



(74) 



Imaginary part of this expression is given in Eq. (13). 



V. EDGE SINGULARITIES 

In the previous section we found that the spectral func- 
tion diverges at the hole mass-shell, which coincides with 
the edge of the spectrum e = ^j, at /c < 0. Here we con- 
sider the behavior of the spectral function in the vicin- 
ity of the remaining kinematic boundaries. As discussed 
in Section II, in all cases the spectral function exhibits 
a power-law suppression at the edge. The suppression 
originates in the phase space constraints which lead to 
vanishing of the spectral function linearly with the dis- 
tance to the edge. Interactions modify the exponent via 
a mechanism analogous to the X-ray edge singularity in 
metals^^: the transition amplitudes in Eqs. (14) and (21) 
acquire a power-law dependence on the distance to the 
edge. In this Section we develop a technique to account 
for this dependence. 

The edge of the spectrum corresponds to final states of 
the transition in which all of the momentum and energy 
are carried by a single hole, see Section II. At energies 
close to, but not precisely at, the kinematic edge, the final 
states may contain, in addition to this "deep hole" , an 
arbitrary number of low-energy particle-hole pairs near 
the two Fermi points. 

When the distance in energy to the corresponding 
edge is small compared to k'^/m, the hole can be treated 
as distinguishable from the rest of the particles in the 
system^'^°. Formally, this amounts to projecting out all 
states except those in the narrow stripes of momenta r, I 
near the Fermi points, and a strip d deep below the Fermi 
level. (Note the difference between the narrow c?-subband 
defined in this Section and the wide one in Section IV). 
The d-subband hosts a single hole, and the interaction 



of this hole with the rest of the system results in the ex- 
citation of the particle-hole pairs in the subbands r and 
I. 

We consider first the behavior of the spectral function 
near the hole edge e « at A; > 0, and then proceed 
with the consideration of the particle edge e w ^/j. 



A. e ^ ^fc , fc > 

In this case the deep hole in the final state of the transi- 
tion is located at momentum close to k (relative to —pp) 
on the left-moving branch of the spectrum, see Fig. 10. 
In order to describe the final-state interaction, we project 
the Hamiltonian (34) onto narrow strips of momenta 
shown in Fig. 10, ffeff = VHV, where V projects onto 
states within r,l, and rf-subbands, while the remaining 
states are regarded as either occiipicd or empty. In or- 
der to extract the dependence of the spectral function 
on the distance to the edge with the leading logarith- 
mic accuracy, it is sufficient to carry out the projection 
to zero order in interactions for the Hamiltonian and to 
the lowest nonvanishing (first) order for the observable. 
The omitted higher order contributions would contribute 
beyond the leading logarithmic accuracy. 




FIG. 10: Solid lines indicate the states included in the ef- 
fective Hamiltonian (76)-(80). The d-subband hosts a single 
hole. The interaction of the hole with r and /-subbands leads 
to the excitation of low-energy particle-hole pairs. Production 
of multiple particle-hole pairs results in a power-law depen- 
dence of the spectral function in a close analogy with X-ray 
edge singularity in metals. 

We introduce slowly varying in space fields 



\k\<ko 



|fe'|<fec 



(75) 

where fco "C fc is the high-momentum cutoff, and write 
the effective Hamiltonian in the coordinate representa- 
tion as 



HeS = Hrl + Hd + -ffint- 



(76) 



Here 



Hri = jdx[r\-ivd^)V +i^^\ivd^)i^^] (77) 
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and 



(78) 



describe r, I, and d subbands, respectively. In writing 
Eqs. (77) and (78) we linearized the spectrum within the 
respective subbands, so that 

^k+k' =^k- Vdk', Vd = v- k/m. (79) 

The third term in the right-hand side of Eq. (76) de- 
scribes the interaction between the subbands, 



j dxp%Vo - Vk)p' + (Fo - y2p,-fc)p1 
+ [Vo - V2p^) j dxp^{x)p\x), (80) 



where p^'^x) = tp^'''^ {x)ip^'''{x) are particle densities in 
r, Z-subbands, and p'^{x) = il)'^{x)'4)'^\x) is density of 
holes in d-subband. In writing Eq. (80) we have set 

K = Uq. 

Our goal is to evaluate the hole contribution to the 
spectral function 



= Re - / 



dte-*(V'f (i)V'^'(O)). (81) 



The projection should now be appUcd to the operator 
ip^ in Eq. (81). In this case, however, the lowest order is 
insufficient, as the state with momentum k on the right- 
moving branch of the spectrum lies outside the subbands 
r,l,d of the effective Hamiltonian. Instead, the relevant 
contribution is generated in the first order in the inter- 
action strength. 

The higher-order contributions can be found with the 
help of a unitary (Schricffer- Wolff) transformation that 
decouples states within the subbands r, I, and d from the 
rest of the system. Consider the following term in the 
original Hamiltonian (34), 



5H 



(82) 



\ki\<ko 



Schrieffer- Wolff transformation H e^^He~^^' elimi- 
nates the off-diagonal contributions such as Eq. (82). To 
the lowest (first) order in the interaction strength, the 
generator of such transformation reads 



Vk-V2p^ 



ki\<ka 



X (l-P)Vf<<V'' 



fe + /C3 



V 



H.C., 



so that [iS, Hq] = —dH, where Hq is the noninteracting 
part of Eq. (34). 



As far as the effective Hamiltonian (76) is concerned, 
the transformation leads merely to the correction to i?int) 

which amounts to negligible second-order corrections to 

the coupling constants in Eq. (80). 

The same Schrieffer- Wolff transformation applied to 
the operator ijj^ yields 



(83) 



Vk - V2p^ 1 



R-tL T. E ^k!i^t^k+ks^k,M+k3- 
ki\<ko 



This contribution, as well as those generated in higher 
orders, conserves separately the numbers of right- and 
left-movers. With this constraint, the first-order contri- 
bution (83) is a product of the least possible number of 
operators acting in the subbands r, I, d, hence, it is the 
most relevant one as far as the behavior of the spectral 
function near the edge is concerned. 

Using now Eqs. (75), (81), and (83), we find 

Ak{e) (xRefdx f dte-''\^''{x,t)^{0,0)) (84) 

J J —CO 



with 



^{x) = ip^^{x) ip^x) il)'^{x). 



(85) 



To proceed further, we bosonize fermions in r and l- 
subbands according to 

^P-'I = v/fc^e±''^"\ [^^'\x),^^\y)] = ±insgn{x - y), 

(86) 

where the upper /lower sign corresponds to r/Z, and ko is 
high-momentum cutoff for bosonic modes, 

(^'■''(x)v5'-'n0)- [¥^'■''(0)]') = -ln(lT27r*fcox). 

In the bosonic representation Eqs. (77) and (80) read 

Hri = — y : idxid^ip'^r (87) 



a=r,l 



/dr 
^ -p" [{Vo - Vk) {Vo - V2p^-k) d.^""] 

+ (^0 - V2p,) j [d^^nM , (88) 

and Eq. (85) becomes 

^{x) = ko e' ['P^(^)+v\^)] (89) 

The correlation function (84) with ^ given by Eq. (89) 

and with the dynamics governed by the effective Hamil- 
tonian (78), (87), and (88) can be evaluated exactly. This 
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W 



can be done by following the steps outlined in Section IV. 
Equivalently, one can diagonalize the effective Hamilto- 
nian by an appropriate unitary transformation^'^'^. The 
generator of such transformation reads 

= \J dxpa{x)[n2pF-kV^{x) - Hkv'ix)] (90) 

+ ^Jdx[^^d,^'-^^d,^^] 

with fXk introduced in Eq. (9) above. In writing Eq. (90) 
we omitted all but the first-order in interactions contri- 
butions (recall that second and higher order corrections 
have been neglected in the derivation of the effective 
Hamiltonian). The transformation decouples the sub- 
bands from each other to linear order in interactions, 

e'^H.se-'^ = Hri + Ha, (91) 
and also modifies the operator ^, see Eq. (89), 

(92) 

Since the transformed Hamiltonian (91) is quadratic, 
evaluation of the correlation function Eq. (84) is straight- 
forward. Keeping only linear in fip terms in the expo- 
nents, we find 

(*t(a;,f)^(0,0)) oc 5{x + Vdt)e^^* (93) 

X fco [1 + 2mko{vt - x)] -1-'^^''--'=+''=''^' 

X fco [1 + 2mko{vt + x)] . 

Substitution of Eq. (93) into Eq. (84) then yields the 
spectral function near the edge. 



Ak{e) oc (^fe - e)^ 



■fJ'k-\-fJ'2pjp-k—'^(J'2pp 



0{^k - e) , (94) 



see Eq. (28). This result is valid when e is close to the 
edge, — e <C k^/m, and for < A; < 2pi? (note that 
for A; > the above derivation should be modified as in 
this case d-subband in the corresponding effective Hamil- 
tonian belongs to the right-moving branch of the spec- 
trum). 

The interaction- induced corrections to -1-1 in the ex- 
ponent in Eq. (94) (which comes from the phase-space 
constraints) allow for a simple interpretation based on 
the analogy with X-ray edge singularity^^ in metals: —fik 
originates in the attractive interaction of the deep hole 
on the left-moving branch of the spectrum with a soft 
left-moving particle; H2pF-k is due to the repulsion be- 
tween the deep hole and a soft right-moving hole; fi- 
nally, —2^2pF is due to the attraction between a soft 
left-moving particle and a soft right-moving hole (note 
that both particles are soft, hence the factor of 2). 

For small k, the exponent in Eq. (94) simplifies to 

1 - Mfe + M2pF-fe - 2/i2pF ^ 1 - Mfe - /"2pF ^ 1 - 270, 

where we used /i2pi? = 270 and and the fact that for a 
generic interaction /x/c ^ when k ^ 0. 



B. e^|fe, fe<0 

The consideration in this case is very similar to that 

in Section VA. The effective Hamiltonian accounting for 
the final-state interaction again consists of subbands r 
and I near the two Fermi points, and rf-subband (now 
centered at momentum — fc > 0) deep below the Fermi 
level on the left branch of the spectrum, see Fig. 11. The 
projection, carried out to the lowest order in interactions, 
yields Eqs. (76)-(80) with the replacement £,k — Cfc ^^id 
Vd-^v + k/m in Eq. (78), and A; ^ -fc in Eq. (80)7 



A;<0 



V 



-Pf 
-PF-k 



Pf P 



FIG. 11: Solid lines indicate the subbands r,l, and d in the 
effective Hamiltonian for the evaluation of the particle contri- 
bution to the spectral function at e ^ ^fc, fc < 0. 

The particle contribution to the spectral function is 
given by 

1 /■OO 

Ak{e)=Be- dte-*<^f (i) (0)). (95) 

Similar to above, application of the Schrieffer- Wolff 
transformation to the operator tp^'' yields the relevant 
contribution in the first order in the interaction strength, 



Vk — V2pp+k 1 



1 ^k^^t^-k+ks^k^+k^M- 



\ki\<ko 



(96) 

When written in the coordinate representation, Eqs. (95) 
and (96) give 

Ak{e) (xReJdxJ^ dte"'\^{x,t)^\0,0)) (97) 



with 



*t(x) = '(/''■■^(a;) V'^(x)V'''(x). 



(98) 



The spectral function (97) is evaluated by bosonizing 
the Hamiltonian and diagonalizing it by a unitary trans- 
formation, just as it is done in Section VA. This proce- 
dure yields 

Ak{e) oc (e-^)i-M'e-M2p^+'=+2M2p^0(e_^) . (99) 

This result is valid at e — ^ k'^/m, and for all k in the 
range —2pF < k < 0. Note that /U2pp and I22pp+k enter 
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the exponent in Eq. (99) with opposite signs compared to 
those in Eq. (94). This is because the interaction of the 
deep hole with a particle on the r-branch is attractive, 
whereas the corresponding correction to the exponent in 
Eq. (94) originates in the repulsion between the deep hole 
and another hole in r-subband. 

As expected, the exponent in Eq. (99) is invariant upon 
the momentum inversion fc <-> 2pi? + fc. At very small |A:|, 
the exponent in Eq. (99) simplifies to 

1 - M/s - M2pF+/s + 2/i2pF ^ 1 + M2pF = 1 + 270. 



C. e-»Cfe, fe>0 

In this limit the final state of the transition involves 

creation of a hole with the momentum — fc on the right- 
moving branch of the spectrum, see Fig. 12 and discus- 
sion in Section II. In order to evaluate the spectral func- 
tion in the leading logarithmic approximation, it is suf- 
ficient to consider the effective Hamiltonian consisting 
of d-subband hosting the hole, and r-subband to allow 
for the creation of low-energy particle-hole pairs near the 
right Fermi point, see Fig. 12. (Inclusion of left-movers 
would merely add a second-order in interaction correc- 
tion to the exponent.) Carrying out the projection and 
linearizing the spectrum within the subbands, we find 



with 



Hr = J dx tp'^^ {—ivdx) tp^ , 
(here Vd = v — k/m), and 



(100) 

(101) 
(102) 



{Vo-Vk) dxp^{x)p^{x). (103) 



fc>0 



i 



r 



Pf- 



Pf P 
-k 



FIG. 12: States included in the effective Hamiltonian (100)- 
(103) for the evaluation of the spectral function at e — > ^fo 
fe > 0. ~ 

After bosonizing Eqs. (101) and (103), the effective 
Hamiltonian can be diagonalized by a unitary transfor- 



mation with generator 

W = -^ J dxp'^ix) ^"^{x). 



(104) 



The transformation yields, to first order in the interaction 
strength, 



(105) 



where Hr = {v / Air) j dx {dx^^Y ^"^^ is given by 
Eq. (102). 

Since the right-moving state with momentum k ^ ko 
lies outside the domain of the effective Hamitlonian (100), 
application of the lowest-order projection to the operator 
fp^'' , see Eq. (95), is insufficient. Similar to above, the 
relevant contribution emerges in the first order in the 
interaction strength and reads 



\hi\<ko ^'^ ^'>-fc 

Since fco <IC fc, we can expand here 



VV'feJV'_fe + fe3 4l+fc2,fe3- 

(106) 



Vk-ki - Vk-k2 ~ (^2 - ki) 



dVk 
dk 



Passing over to the coordinate representation, we find 
that the spectral function is given by Eq. (97) with 

*+(x) = ip''^{x){-idxW^ix) i)^{x) . (107) 

Bosonizing Eq. (107) according to Eq. (86) and taking 
proper care of the point splitting (see Appendix C for 
the details), we find 



*t(a;) = 27rfc^e-2''^''(^)V'^(a;). 



(108) 



The unitary transformation (104)-(105) that diagonalizes 
the effective Hamiltonian also modifies the operator 

e^^^ix)e-'^ (X e-'(2-'''=/2)^''(^V^(a;). (109) 

With the dynamics governed by the quadratic Hamilto- 
nian (105), it is straightforward to calculate 

(* (x,t)\E'^(0,0)> cx Six - Vdt) e-'^* (110) 
xk^[l + 2Triko{vt-x)] 

Substitution of Eq. (110) into Eq. (97) then yields 

Ak{e)cx{e-^f-'"'d{e-^) (111) 

for the spectral function at e — ^fe <C k-^/m. The fact 
that the interaction-induced correction to the exponent 

is — 2/ife is due to the presence of two soft particles inter- 
acting with the deep hole, contributing —pk each. 
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VI. CALOGERO-SUTHERLAND MODEL 

In this section we consider a solvable model of interact- 
ing ID fermions, the Calogero-Sutherland (CS) model. 
The corresponding Hamiltonian in the first-quantized 
form reads^* 



H 



y- 

^ O r. 



2m dxj 



Y^V{x,-x,), (112) 



where V{x) is a periodic version of the inverse-square 
interaction potential, 



V{x) 



\{\~l)/m 
{L/n)^sm'^{nx/L) 



(113) 



Correlation functions of the CS model exhibit rather 
unusual behavior. For example, the dynamic structure 
factor S{q,uj) differs from zero in a finite interval of fre- 
quencies®'"'^^'^^'^^, just as it is for free fermions. On the 
contrary, for a generic interaction the structure factor 
has a high-frequency "tail" 5* oc q'^/oj'^, which emerges 
already in the second order of the perturbation theory in 
the interaction strength^'^. However, close examination 
of the corresponding perturbative formula for the struc- 
ture factor (see Eq. (18) of Ref. [6]) indeed shows that it 
yields zero for Vk cc \k\. 

Similarly, substitution of C/fe = Vk oc |fc| into Eqs. (40)- 

(43) yields l/2rfc = -ImE^'*^(^fe) = 0. Just as it is 
the case with the absence of the high-frequency tail in 
S{q^Lo), the apparent vanishing of the relaxation rate in 
perturbation theory suggests that 1/rfc = is the ex- 
act relation for the CS model. Indeed, this agrees with 
the exact results for the Green function^^ obtained for 
specific values of A. 

Vanishing of 1 /t/j for particles is a peculiar property of 
the CS model related to its integr ability^ ^. Therefore, in 
the context of this work, we concentrate on the hole con- 
tribution to the spectral function; we believe that the CS 
model results for the hole (e < 0) region of the spectrum 
are generic. 

Single-particle correlation functions for the CS model 
have been studied extensively^^'^° in the context of the 
exclusion statistics^^ . Such interpretation is possible be- 
cause the inverse-square potential is impenetrable. Re- 
quiring the many-body wave function to obey a certain 
symmetry with respect to the permutations of the parti- 
cles' coordinates amounts to merely choosing a rule ac- 
cording to which the wave function is assigned an overall 
phase. For CS model, the phase depends on the order- 
ing of particles, but not on their coordinates. Operators 
which do not permute particles (e.g., the local density 
operator) do not affect the phase factor. Accordingly, 
the statistics of bare particles is immaterial as far as the 
evaluation of, say, the dynamic structure factor S{q,uj) 
is concerned^". 

However, the situation with single-particle correlation 
functions is more subtle. The anyon creation and annihi- 
lation operators constructed and studied in Ref. [29,30] 



describe fermions only for odd integer values of A. Rather 
than attempting to derive fermionic Green function for 
general values of A, below we will use the results of 
Refs. [30] to evaluate Ak{e) for e < and A — odd integer, 
and then employ an analytical continuation to extend the 
result to arbitrary values of A. 

The excitations of CS model can be described in terms 
of quasiparticles and quasiholes. Quasiholes are charac- 
terized by fractional inertial mass fh = m/X and veloci- 
ties Vi in the range \vi\ < v = pp/rh, where pp — Trn is 
the Fermi momentum {n is particle concentration). On 
the contrary, quasiparticles have velocities \vi\ > v and 
their mass coincides with the bare mass m that enters 
Eq. (112). 

Consider now a state obtained by a removal of a single 
particle (with mass m) from the ground state |0). For 
integer A, this change of mass can be accommodated by 
creation of exactly A quasiholes. It turns out that this 
simplest possibility is, in fact, exhaustive'^''. 




FIG. 13: Support of the spectral function in (p, e < 0) half- 
plane. The upper boundary of the region where A ^ com- 
prises of A identical parabolic segments (A = 5 in the figure). 
We concentrate on the low-energy sector k — p — pp ^ (see 
the inset). 

The region of support of Ap{e) in {p, e)— plane can now 
be deduced from the energy and momentum conserva- 
tion. Indeed, using the well-known expressions^^ for the 
energy and momentum of a state with A quasiholes with 
velocities Vi, one can write 



A _ 

Em 
y 

1=1 



A 



{v^-vf) = 0, p-Vmv, = 0. (114) 



For odd integer A, these equations have a solution for 
{vi} provided that e and p lie within the shaded region 
in Fig. 13. The upper boundary of this region, the solid 
line e — ^p, comprises of A identical parabolic segments, 

^=^[iP~^lPFf -pI], \p~21pf\<pf (115) 

with integer Z, \l\ < (A— 1)/2. The support of (e) is also 
bounded from below by the dashed line e = p^/2m — ep 
with ep = mv'^/2; this boundary would be absent for 
generic values of A. 

We concentrate here on the low-energy sector with p w 
Pp, \e\ <C £p (see Fig. 13) and "shift" the momentum 
according to k = p — pp. The dependence of Ak{e) on e 
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has a threshold, 

Ak{e)(xe{^-e), a = -u|fc| + fc72m, (116) 

and our goal here is to find the behavior of A at a fixed 

k when e approaches the threshold. This can be done by 
writing the spectral function in the form of a multiple 
integral, 



/y A 
Y[dviF{{vi}) 



(117) 



xS[k+pF-mY: Vi] S[e+ (m/2)E(«' " ■ 

The (5— functions here reflect the conservation of mo- 
mentum and energy, see Eq. (114), and the form-factor 

F oc\{{vi}\tpk\0)f was found in Ref. [30], 



F{m = n 



nc 



2 2\ 
V — Vi ) 



(118) 



1<«2 



Eq. (117) is nothing but the Lehmann representation of 
the spectral function, with the final state of the transition 
\{vi}) parametrized by the velocities of A quasiholes. 

For <^ pf and ^fe — e ^ ep, the final states {vi} 
contributing to A have n\ = {\ — l)/2 quasiholes with 
velocities Vi ~ —v and nx + 1 quasiholes with velocities 
Vi « +v. It is therefore convenient to introduce new 
variables 



Vi = 



V + Vi 



V 



After expansion of the form-factor to the lowest nonvan- 
ishing order in Xi,yi <^ 1, Eq. (117) takes the form 



Ak{e) oc — TT/ dxi dyj fn^+i{{xi}) fn^{{yj}) 
ep -^tJo Jo 

(119) 



X S 



X S 



^+^1-1(^2+12) 

2mv^ 4 
e — vk 



+ Yi--{X2+Y2] 



Here 



"A-l-l 



i=l 



and /jv({-2i}) is a function of N arguments Zi, . 
given by 

ii<i2 i 

this is a homogeneous function of degree 

l{[2(iV-n.)-l]^-A^}. 



(120) 



Cat 



As it is easy to check, 

X2 ~ Xf < Xi < 1, Fa ~ « Fi < 1. 

Moreover, when e is relatively far from the threshold, 
at — e 3> k'^/2fh, Xi and Yi are of the same order, 
Xi ^ Yi. In this limit X2 and I2 in the arguments of 
the (5-functions in Eq. (119) can be safely neglected, after 
which the integration reduces to a power counting which 
yields 



Ak{e) oc \e + vky° \e - vk\ 



|7o-l 



with 



4A 



a-e»A:7m (121) 



(122) 



Eqs. (121) and (122) reproduce the standard Luttinger 
liquid result. At \vk — e| <C v\k\ Eq. (121) agrees with 
Eq. (4) above. Since the exponent 70 is an analytical 
function of A, Eq. (122) is valid for all A rather than for 
integer values only. Indeed, in the weak interaction limit 
|A — 1| <C 1, Eq. (5) yields 70 ~ (A — l)/2, in agreement 
with the corresponding limit of Eq. (122). 

The behavior of the spectral function in the immediate 
vicinity of the threshold, at ^fe — e < k^/fh, requires more 
delicate consideration. Indeed, in this limit the velocity 
of one of the quasiholes approaches 



vo{k) = — ■ysgn(fc) -|- k/m, 



(123) 



while all other quasiholes have velocities — > iv. In 
other words, close to the threshold, almost all of the mo- 
mentum and energy are carried by a single quasihole, 



m 



In order to evaluate the integral in Eq. (119), we note that 
for ^fe — e <C k'^/fh and a = \k/pp\ <C 1 the conservation 
laws [i.e., the (5— functions in Eq. (119)] imply that for 
k<0 



X, 



Actually, for A: < it is just one of x,, say, a;„^+i, which 
is close to a, while for i = 1,. . .n\ one has ^ a^. It is 
therefore convenient to write 

and introduce Xi = Xi — a = ~ a^. In the 

leading (second) order in a the ^— functions in Eq. (119) 
can be approximated by 



6 



2mv'^ 



X, 



Yi 



[note that < (Cfc ^ e)/mv'^ < a^]. At the same time, 
Xnx+i in the form-factor should be replaced by a. With 
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these approximations, the remaining integrations are eas- 
ily carried out resulting in 

^fc<o(e)cx(a-e)"'°"\ 0<a-e<fc7TO (124) 

The fractional part of the exponent here is twice that 
in the Luttinger liquid limit, see Eq. (121). This is in 
agreement with the leading logarithmic approximation 
result Afe<o(e) oc {^k — e)'''""^, see Eq. (13). Indeed, sub- 
stitution of Vk = — A(A — l)Tr\k\/m into Eq. (9) gives 
Mfe = A*o ~ A — 1, independently of k (this agrees with 
the exact value^° of the exponent that governs the diver- 
gence of the structure factor, /xq = 1 — 1/A). Eq. (8) then 
yields 7^ = Mo/-^ ~ 27q. Note also that the position of 
the spectral edge £,k, see Eq. (116), in the limit A ^ 1 
agrees with Eq. (24). 

Evaluation of Ak{e) for fc > proceeds similarly and 
yields 



^fe>o(e) oc (^k - e) 



< a - e < fc'M (125) 



with 71 given by 



2 (A -3)2 
7i = 



4A 



(126) 



For a weak interaction 7o + 7i ~ 1 — 2/io, in agreement 
with Eqs. (28) and (94). Note that the exact exponent 
in Eq. (125) is positive for any strength of interactions. 

We have thus demonstrated that, at least in the hole 
region of the spectrum e < 0, the behavior expected for 
a generic ID system with a nonlinear dispersion is con- 
sistent with the exact results obtained for the Calogero- 
Sutherland model. It would be interesting to find a solv- 
able model for which our conclusions for the particle re- 
gion of the spectrum can be similarly tested. 



VII. CONCLUSIONS 

The total number of particles and the total momen- 
tum are good quantum numbers for an isolated homo- 
geneous fermionic system, regardless of its dimensional- 
ity and the interaction strength. In higher dimensions 
{D > 1) and for moderately strong interactions, an ex- 
cited state of the system with one extra particle and with 
momentum p is rather similar to the corresponding state 
of a free Fermi gas. The similarity is encoded in the en- 
ergy and momentum dependence of the spectral function 
Ap{e) = — (I/tt) ImGp(e). The spectral function satisfies 
the exact sum rule, 



j deAp{e) 



1, 



(127) 



In a Fermi liquid [D > 1), the sum rule is almost com- 
pletely exhausted by the Lorentzian Eq. (1). The cor- 
responding peak is centered at the quasiparticle energy 
and has the width 1 /2rp which decreases with the in- 
crease of the Fermi energy, see Eq. (2). On the contrary, 



in a Luttinger liquid the spectral function, see Eq. (4), is 
manifestly non-Lorentzian. It diverges on the mass-shell 
e = ^fc, vanishes at e < ^fc, and decays slowly with e at 
e > 6- 

We demonstrated that for a nonlinear dispersion re- 
lation with positive curvature, see Eq. (6), the domain 
where Ak.{e) differs from zero extends below the mass- 
shell, e > ^fc — k^/m. The mass-shell e = ik now falls 
within the domain of continuous spectrum. In this situa- 
tion the Luttinger liquid resTilt Eq. (4) is no longer valid. 
Instead of a one-sided power- law singularity at e — ^ ^fc, 
[cf. Eq. (4)], the spectral function has a Lorentzian peak 
centered at the mass-shell. The width of the peak l/r^, 
see Eq. (10), is much smaller than that in higher dimen- 
sions. In one dimension, the finite relaxation rate l/r^ 
emerges only in the fourth order of the perturbation the- 
ory in the interaction strength and increases with k as k^ 
(here k is momentum relative to pp). At large enough 
momenta k, which are still exponentially small in the 
inverse interaction strength 



k> ppe 



(128) 



the Lorentzian peak in the spectral function carries most 
of the spectral weight. In other words, for fc satisfying 
the inequality (128), the dispersion nonlinearity restores 
the Fermi liquid character of the particle part of the spec- 
trum. The numerical coefficient a in the criterion (128) 
depends on the fraction of the spectral weight chosen to 
fall within the Lorentzian. 

It should be emphasized that dispersion nonlinearity 
brings about particle-hole asymmetry into the problem. 
Indeed, we found that the nonlinearity does not affect 
qualitatively the hole region of the spectrum e < 0. The 
spectral function here resembles that in the Luttinger 
liquid, although with the modified value of the exponent. 

Finally, we mention that the emergence of a finite 

quasiparticle lifetime brings some ramifications for the 
theory^ of the structure factor S{q,uj) of an interacting 
ID fermion system. It smears the nonanalytical behav- 
ior of S{q,Lo) at to = vq + q'^/2m. The smearing, how- 
ever, affects only a tiny portion of the frequency interval 
I'*-' ~ ^ q^/2ni. Indeed, the width of the smearing 
region scales with q as I/tq oc q^/rn?, while the width of 
the interval of interest is g /2m. 
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APPENDIX A: VERTEX CORRECTIONS 

In this appendix we consider the dimensionless cou- 
pling vertex r characterizing the interaction of the d- 
particle with the bosonic field </> describing excitations 
within the subbands r and I, see Eq. (64). Our con- 
cern here is whether corrections to T due to interactions 
within the d-subband affect the leading logarithmic se- 
ries considered in Section IV. The vertex depends on 
four variables, F = r(e, k, oj, q). The dependence of F on 
its arguments around the mass shell is not singular. It is 
clear then that the relevant object is the on-shell vertex 
correction 

6T = T{^k,k,±vq,q)-l 

(the bare value of the vertex is 1). 

By construction of r, Z-subbands (see Section IV A), 
the momentum q transferred to the bosonic field <p is re- 
stricted to \q\ < Xk. Accordingly, the relevant frequencies 
CO = ±vq satisfy < Xvk <C e. It then follows that the 
two d-particle Green functions adjacent to the bosonic 
field (j) in the diagram such as that shown in Fig. 9(b) 
carry frequencies of the same sign, sgn(e) = sgn(e + oj). 
The on-shell vertex correction can then be expressed via 
the time-ordered d-particle self-energy (e) as 



Eq. (72). The coordinate integration can be performed 
at once. With the help of the convolution theorem, the 
remaining integral over time be written as 



mv 



rlf' 

— Gt{e-e')Tk{e'), (Bl) 



where 



(B2) 



is the bare d-particle Green function, see Eq. (59), and 
^fe(e) is given by 



.Ffe(e) = /dte"* [l+i (yd + v)t 

mv 



X (1 - iXk{vd - v)t) 



-M=_/4 



(B3) 



where we denoted fi- = /i^ and = ii2pp+k for brevity. 
The integration in Eq. (B3) is carried out using the con- 
volution theorem and the relation 



5T 



de 



+ 0{q). 



(Al) 



The latter is discussed in details in Sections III and 
IV above. In the second order m. d — d interaction 
the self-energy is purely real and is given by S^^^ <^ 
7o(e - ^fe)ln(A/A), see Section IV. Eq. (Al) then yields 



ReJF ~ 7^1n(A/A) «; 1, 



(A2) 



where we used Eq. (57). In the fourth order the self- 
energy acquires a finite imaginary part on the mass-shell, 
ImE^'*''(^fc) = — l/2Tfe. Taking into account that ImE^^' 
varies with e on the scale 
we find 



k /m, and using Eq. (Al), 



Im^F ■ 



l/Tk 

k'^/m 



< 1. 



(A3) 



According to Eqs. (A2) and (A3), vertex corrections 
due to interactions within d-subband are small. As far as 
the summation of the leading logarithmic contributions 
in Section IV is concerned, all such corrections can be 
safely neglected. 



APPENDIX B: DERIVATION OF EQ. (72) 

In this appendix we supply technical details needed 
to perform the Fourier transform of Eq. (71) leading to 



/ 



dte 



27r e{±e) 



(B4) 



where k ^ Xk'^ /m. Taking into account that vj, — v 
k/m and v^ + v '^2v, we obtain for /ij- <C 1 



e'(e' + e) V 



(B5) 



Eq. (B5) is valid for \e\ < k = Xk'^ /m. In the case e > 
we have 



,2 ,,2 



8e 



(1) 'B{,.l/A,l-,l) 



(B6) 



where = see Eq. (8), and B is the standard 

beta-fimction. For weak interactions, ii± <^ 1, one finds 
B « 4:/ lit- The case e < is analyzed in a similar way. 
As a result, we obtain 



me) = 



2e 



7fc 



[fJ,le{e) + ^J,^_e{-e)], (B7) 



covering both cases. We now substitute Eqs. (B2) and 
(B7) into Eq. (Bl). Rewriting the integral in the latter 
equation as the contour integral we find 
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/ \p 

Gk e;—,Ak 
\ mv 



-Ik 



dz 



(-^) 



47fc Jc+ 27ri e - ^fe + i/2rfe - ^ 47^ 2m e - ^fe + V2Tfe - 2: 



dz 



(B8) 



where the function is chosen to have a cut in the 
complex plain along the negative real axis and contours 
C± run clockwise around the positive (negative) half of 
the real axis. Closing the contours at infinity and taking 
the residue at 2: = e — ^/j + i/2rfe we arrive at Eq. (72). 



APPENDIX C: DERIVATION OF EQ. (108) 

Consider the operator 

^\x) = ^f\-id^Y\r' (ci) 

entering Eq. (107). Proceeding in a standard fashion, we 
subsitute here '>lf = ^/koe'^'^ with (p = (p'^, see Eq. (86), 
and write the derivative as a finite difference. 



-ikt 



2A 



■0 ^-iip(x) 



^-i(p{x+A) _ ^-itp(x-A) 



(C2) 



On making use of the identity 



gag5 ^ . ga+5. g(a6+i (a=+6^)) 



I 

valid for any operators a and b which are linear in bosonic 
fields (the colons denote the normal ordering), we find 

X (l±27rzA:oA). 

We now approximate here 

-i[p{x)+ip(x±A)\ . g-2i<^(a:). ^ ^2{v'^ (0)) - 2iip(x) ^ 

SO that 

g-iy(x)g-i¥'(x±A) _ g-2i<^(x) _^ 27rifcoA) . (C3) 

Eqs. (C2) and (C3) yield 

$t(x) =27rfc^e-2'^(^\ (C4) 

Substituting this expression into \E'''' = '^^ijj'^ [see Eqs. 
(107) and (CI)], we arrive at Eq. (108). 
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